AD.4  037694 


AFCEC-TR-76-47 


PREDICTION  OF  TEST  CELL  VISIBLE  EMISSIONS 


DIRECTORATE  OF  ENVIRONICS 


D D C 

n<  4 1977 


Ik 


DECEMBER  1976 


FINAL  REPORT  : JUNE  1976  - NOVEMBER  1976 


E 

A 


Approved  for  Public  Release;  Distribution  Unlimited. 


AIR  FORCE  CIVIL  ENGINEERING  CENTER 

(AIR  FORCE  SYSTEMS  COMMAND) 

TYNDALL  AIR  FORCE  BASE 
FLORIDA  32401 


JfJCLAbb  l Fib  LI 


Cl  I Ml  T Y CL  ASCI  Kit  ATION  0>  ' HlS  PAUL  (H'h«ii  Dm  a l-.iHoretll 


REPORT  D(  CUMENTATION  PAGE 


f i~  nrAn»t  miwnr 

1 AFCEC-'.Vtt-  76 


CK5SI0N  NO 


TITLE  (nrui  Subtitle) 

(/  7 ^PREDICTION  OF  TEST  CELL  VISIBLE  EMISSION^ 


Samuel  ,.P.  /Tinch 


/ I n i in 

( Ub 

Arland  M . 


Major,  USAF 
CO) 

Cap':  , USAF,  HSC 


JL.  PT  |l  p nUKIili  fiin-frTtTT  ATrorTTAMF  and  address 

Air  Force  Civil  Engineering  Center  k 

Environics  Directorate 

Tyndall  Air  Force  Base  FL  32403 


IL  CONTROLLING  OFFICE  NAME  AND  ADDRESS  , . I ’U 

Air  Force  Civil  Engineering  Center  ( ’ fh 
Tyndall  Air  Force  Base  FL  32403  V — 


.aOl&V-NAME  A AOORESSf/I  different  from  ContrnIHnt  Office) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


NT'S  CATALOO  NUMBER 


■ CONTRACT  OR  GRANT  NUMBERf*; 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  «,  WORK  UNIT  NUMBERS 

63723F 
2103  7A29 


7 TiW™^ mop i u a rf” 

Dec  £8  76 


1 


13  NUMBER  PAGES 

41 


IS.  SECURITY  CLASS,  fo ( thl,  report) 

UNCLASSIFIED 


IG.  DISTRIBUTION  STATEMENT  (ol  (Ala  Report) 


15a.  OECLASSlFlC  ATION/ DOWN  GRADING 
SCHEDULE 


Approved  fnr  pnb.1- 


tribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  <h»  sb»trmct  ent«r»ci  In  Block  20,  1 1 dl  H*r«nt  from  Raport) 


18.  SUPPLEMENTARY  NOTES 

Available  in  DDC . 


19.  KEY  WORDS  ( Contlnua  on  roverao  *ide  It  nacAAAary  And  Identity  by  block  numbar) 


Environmental  effects 
Jet  engine 
Air  pollution 


Plumes 

Cell,  emission 
Carbon  particles,  haze 
Air  flow  rates,  cell 


C.  ABSTRACT  (ConUnu*  on  t over  a a a/de  It  ne  coanary  And  IdantUy  by  block  numbar) 

A theoretical  correlation  between  jet  engine  test  cell  plume 
opacity  and  Society  of  Automotive  Engineers  (SAE)  engine 
smoke  number  (SN)  was  developed.  Intermediate  results  provide 
soot  loading  at  the  exhaust  plane  of  the  engine,  engine  exhaust 
flow  rates  and  test  cell  total  air  flow  rate.  The  specific 
particle  extinction  coefficient  which  correlates  light  scattering 
properties  with  soot  loading  is  the  most  difficult  parameter  to  _ 


DD  i j an  *7i  1473 


EDITION  OF  1 NOV  66  IS  OBSOLETE 


UNCLASSIFIED 


foo^e 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wim  Dele  Holered) 


390  9Zdr 


'iA 


UNCLASSIFIED  

SECURITY  CLASSIKICATION  OE  THIS  PAOEpthwi  Dtls  Bnltrtd) 


fO 


nt 


20  (continued) . 

define.  A value  predicted  from  Mie  theory  is  used  for 
correlations  close  to  the  exhaust  plane  of  the  engine.  At 
the  exit  plane  of  the  exhaust  stack  agglomeration  and  scouring 
change  the  particle  size  distribution  and  individual  particle 
density  so  that  theoretical  prediction  is  difficult.  The 
value  for  specific  particle  extinction  coefficient  at  the 
exhaust  stack  was  chosen  based  on  what  little  empirical  data 
was  available.  When  additional  data  becomes  available,  it 
should  be  possible  to  define  this  parameter  more  precisely. 
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SECTION  I 


ESTIMATION  OF  TEST  CELL  PLUME  OPACITY 


To  insure  Air  Force  facility  compliance  with  applicable 
federal,  state,  and  local  air  pollution  regulations,  it  is 
necessary  to  estimate  the  opacity  of  plumes  emitted  by  jet 
engine  test  cells.  Since  Ringelmann  readings  are  subjective 
and  somewhat  inaccurate  (Reference  1) , and  the  costs  of  in- 
stalling  a transmissometer  or  making  visual  readings  on  each 
facility  is  high,  a theoretical  method  using  known  engine 
and  test  cell  parameters  is  desirable.  There  is  no  known 
work  which  relates  engine  particulate  emission  taken  at  the 
exhaust  plane  of  the  engine  to  the  plume  opacity  at  the  exit 
plane  of  the  test  cell  exhaust  stack.  However,  there  are 
several  works  that  deal  with  intermediate  portions  of  the 
required  calculations. 

The  smoking  tendency  of  jet  engine  exhausts  is  most 
commonly  reported  in  terms  of  smoke  number  (SN)  computed  in 
accordance  with  the  Society  of  Automotive  Engineers  ARP  1179 
(Reference  2) . Several  published  works  have  related  the  SN 
to  the  optical  properties  at  the  exhaust  plane  of  the 
engine  (References  3 to  7) . These  have  resulted  in  good 
correlations  with  the  experimental  results  reported  by 
Shaf fernocker  and  Stanforth  (Reference  5) . Wood  (Reference 
3)  has  done  the  latest  work  in  this  area,  and  his  work  is 
the  basis  for  present  relationships  correlating  the  soot 
loading  at  the  exit  plane  of  the  engine  with  opacity  using 
the  Beer-Lambert  Law. 

In  each  test  cell  there  is  substantial  additional  air 
flow  around  the  engine  that  dilutes  the  emissions.  To 
predict  the  opacity  of  the  plume  as  it  leaves  the  exhaust 
stack  of  the  test  cell,  it  is  necessary  to  know  this  air 
flow.  The  air  flow  can  be  calculated  if  the  cell  geometry 
and  the  pressure  loss  caused  by  the  air  flow  through  the 
cell  is  known. 

As  the  particulate  emissions  move  from  the  engine  to  the 
exhaust  plane  of  the  cell,  they  are  affected  by  agglomera- 
tion and  dilution.  Also  at  higher  engine  exhaust  flow  rates 
larger  particles  are  scoured  from  the  walls  of  the  cell. 

This  causes  a change  in  the  mean  particle  size,  the  particle 
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size  distribution,  and  the  density  of  the  individual  particle 
The  light  scattering  properties  of  the  aerosol  are  predicted 
by  Mie  theory  using  particle  size  distribution  and  density. 
Ensor  and  Pilot  (Reference  8)  used  particle  size  and  distri- 
bution to  derive  a K factor  to  facilitate  making  Mie  scatter- 
ing theory  predictions  of  plume  opacity. 

The  purpose  of  this  report  is  to  describe  a method  for 
estimating  test  cell  plume  opacities  based  on  available 
engine  parameters,  test  cell  geometry,  and  test  cell  pres- 
sure depression. 
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SECTION  II 


METHODS 


Correlation  of  plume  opacity  with  test  cell  and  engine 
parameters  involves  four  distinct  elements.  The  first 
three  elements,  (1)  conversion  of  SN  to  soot  loading, 

(2)  applying  engine  combustion  parameters  to  obtain  engine 
flow  rates,  and  (3)  using  test  cell  geometry  to  determine 
total  flow  rates,  provide  the  input  variables  for  use 
in  (4)  the  Beer-Lambert  Law  to  make  the  final  opacity 
estimation.  Each  element  will  be  discusseo  independent- 
ly then  combined  with  the  others  to  make  the  final 
calculation. 

1.  CONVERSION  OF  SN  TO  SOOT  LOADING 


SAE  SN  (Reference  2)  is  currently  the  most  common 
method  of  reporting  the  concentration  of  carbon  particles 
in  the  exhaust  of  a jet  engine.  The  technique  compares 
the  reflectance  of  clean  filter  paper  to  the  reflectance 
of  a filter  that  has  been  stained  by  drawing  a specific 
weight  of  engine  exhaust  gas  through  it.  The  loss  of  re- 
flectance is  a function  of  the  amount  of  soot  deposited 
on  the  paper.  There  have  been  numerous  efforts  to  con- 
vert SNs  to  a soot  loading  in  g/cm3  at  the  exhaust  plane 
of  the  engine,  the  latest  by  Wood  in  1975.  Wood  combined 
the  results  of  previous  efforts  to  derive  an  equation 
that  gives  a good  correlation  with  available  empirical 
data.  His  equation  is j 


SAE  SN 
100 


exp 


h ( 


W 


(1) 


where 

W - weight  of  gas  sample  used  for  smoke 
number  determination  (g) 

A - area  of  filter  spot  (cm2) 

We  - soot  loading  at  the  exit  of  the  engine  (g/cm3) 

Pe  - density  of  the  exhaust  gas  (g/cm3) 

Sp  - specific  projected  particle  extinction  area  (cm2/g) 
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The  constants  b and  C are  empirically  derived. 


For 

SN  £ 

15, 

b = C = 1, 

and 

for 

SN  > 

40, 

00 

o 

II 

A 

and  C 

W/A  is  the  weight  of  the  air  sample  in  pounds  per 
square  inch  of  filter  spot  area.  The  SN  is  normally 
taken  at  a W/A  of  0.0230  lb/in2 (1. 62g/cm? ) . If  no  W/A  is 
reported  in  the  data  for  the  engine,  the  above  value 
should  be  used. 


Sp  is  a number  that  relates  the  reduced  reflectance 
of  the  stained  filter  paper  spot  to  the  mass  loading  on 
that  spot.  There  are  problems  in  finding  a value  for 
this  parameter.  There  is  some  question  as  to  whether  the 
particle  absorbs  light  on  the  filter  paper  the  same  way 
it  would  if  it  were  an  aerosol.  The  staining  of  the 
filter  is  also  known  to  be  nonlinear.  The  error  assoc- 
iated with  this  problem  has  not  been  quantified  but 
appears  to  be  small.  Wood  (Reference  3)  calculated  a 
value  of  8.64  x 101*  cm2  /g  using  Shaf  fernocker ' s and 
Stanforth's  (Reference  5)  data.  A theoretical  value  for 
Sp  has  been  derived  from  the  K parameter  of  ^nsor  and 
Pilot.  Both  Wood  and  Champagne  (Reference  4)  agree  that 
the  value  derived  from  Mie  theory  (6  x 10"cm2/g)  gives 
the  best  overall  results.  This  value  will  be  used  for 
converting  SN  to  soot  loading. 

Wood's  equation  was  rearranged  in  the  form: 


- In  ( 1 


y/ 

SAE  SN.  b C 
100  ' ‘ S 


and 


W 

A 


in  grams  per  cm  was  plotted  against  SN. 


Figure  1 shows  the  results  of  this  plot  using  S_  = 6 x 101* 
and  the  appropriate  values  for  b and  C.  Figure*2 
shows  this  same  plot  with  a straight  line  between  SN  70 
and  SN  15.  The  figure  shows  that  a single  straight  line 
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equation  be  used  to  predict 

(jj  ; 

--  . ~ in  g/cm  for  SN  ■ 15  with  little  error. 

f)e  A 

The  equation  of  this  straight  line  is: 

“s.  W , 10(0.0263(SN)-5.96) 

p A ^ 

e 

Using  the  known  values  of  W/A,  the  ratio,  we/pe,  can  be 
calculated.  This  value,  the  particulate  mass  fraction  in 
the  engine  exhaust,  can  then  be  used  in  conjunction  with 
the  engine  flow  rates  calculated  in  the  next  section  to 
determine  the  rate  of  soot  production  in  grams/sec.  It 
is  felt,  that  engines  producing  SNs  less  than  15  will  not 
excetd  plume  visibility  standards  in  most  test  cells. 
However,  if  a soot  loading  is  required  for  a SN  less  than 
15,  Wood's  equation  (1)  using  b = 1 and  C = 1 will  give 
the  desired  value. 

2.  ENGINE  COMBUSTION  PARAMETERS 

Using  SN  to  derive  a soot  loading  at  the  exit  plane 
of  the  engine  results  in  grams  of  carbon  per  cubic  meter 
of  exhaust  flow.  The  engine  fuel  flow,  fuel  to  air 
ratio,  and  the  stoichiometry  of  the  fuel  combustion 
equations  are  used  to  calculate  the  soot  flow  rate  in 
grams  per  second. 

The  Wright-Patterson  Air  Force  Base  fuels  lab  reports 
JP-4  jet  fuel  is  a mixture  of  8 and  9 carbon  saturated 
hydrocarbons  with  an  average  molecular  weight  of  120 
A.M.U.  The  average  molecular  weight  is  the  result  of  a 
weighted  average  mixture  containing  57  percent  C~  chains 
and  43  percent  Cg  chains. 

The  combustion  of  these  two  chains  must  be  considered 
separately.  Assuming  complete  combustion,  the  equations 
are: 

CgHigU)  + 121*  o2(g)—  8 C02(g)  + 9 H20  (g) 

C9H20(1)  + 14  °2(9)—  9 C02(g)  + 10  H20  (G) 
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The  combustion  of  the  Cg  chain  consumes  12.5  moles  of 
gaseous  oxygen  and  produces6 17  moles  of  gaseous  products. 
This  is  a net  increase  of  4.5  moles  of  products  per  mole 
Cg  fuel  burned.  The  Cg  combustion  results  in  a 5 mole 
increase  per  mole  of  Cg  fuel  burned.  Using  the  weighted 
average  of  (.  57)  (4.5)  + (.  43)  (5.0)  results  in  the  system 
producing  4.7  moles  of  product  gases  per  mole  of  fuel 
burned.  Using  fuel  flow  in  pounds  per  hour,  and  a known 
fuel/air  ratio  (use  total  fuel/air  ratio  for  a turbofan 
engine)  and  a molecular  weight  of  28.95  AMU  for  air,  the 
air  flow  and  the  fuel  flow  in  moles  per  second  are  com- 
puted. The  moles  of  fuel  per  second  is  multiplied  by  the 
stoichiometric  relationship  of  4.7  moles  of  gaseous 
products  per  mole  of  fuel  burned  to  get  the  increase  in 
gaseous  products.  This  is  added  to  the  molar  air  flow  to 
find  the  total  number  of  moles  per  second  flowing  through 
the  engine.  Using  the  average  molecular  weight  of  28.95 
g/mole  for  air,  the  total  exhaust  flow  in  grams  per 
second  is  found.  If  this  is  multiplied  by  the  particulate 
mass  fraction  in  the  exhaust  found  in  the  previous  section, 
the  soot  flow  rate  in  g/m3  is  obtained. 

The  method  described  above  was  used  to  compute  the 
soot  flow  rate  in  the  sample  calculations.  The  result 
was  1.45  g/sec,  and  for  a similar  engine  Grems  (Reference 
9)  reported  1.29  g/sec.  This  is  a fairly  good  agreement 
because  of  relating  particle  extinction  coefficients  on 
filter  paper  to  aerosol  extinction  coefficients. 

3.  TEST  CELL  AIR  FLOW  RATE 

The  total  amount  of  air  flowing  through  a jet  engine 
test  cell  can  be  calculated  by  knowing  the  geometry  of  the 
cell,  and  the  pressure  loss  between  the  air  inlet  and  the 
engine  air  intake  caused  by  the  air's  movement  through 
the  cell.  A typical  test  cell  is  shown  in  Figure  3. 

Three  areas  are  referenced,  and  an  average  velocity 
is  associated  with  each.  Area  B is  the  entry  baffles, 
area  E is  a square  or  mitered  elbow,  and  area  R is  the 
reference  area  where  the  cell  depression  probe  is  located. 

The  A3HRAE  Guide  and  Data  Book  (Reference  10)  shows 
that  incompressible  flow  may  be  assumed  with  less  than  3 
percent  error  when  the  ratio  is  0.95  or  greater. 

Since  reported  cell  depression  pressure  ratios  are 
generally  0.98  or  greater,  the  basic  Bernoulli  equation 
for  incompressible  flow  can  be  used  to  determine  the 
velocity  of  the  air  flow  at  the  inlet  of  the  engine. 
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The  equation  used  is  : 


— 2 P 

VR  S2 

- — + Z_  + + losses 

2gc  R p 


where 

V - average  velocity  in  section  being  considered. 

Z - elevation  above  datum  plane. 

p - density  of  air  at  reported  cell  conditions 
(Reference  11). 

g - is  the  gravitational  dimension  conversion 
*'  constant. 


P - uncorrected  barometric  pressure. 

S1 

P - uncorrected  pressure  at  engine  inlet. 
s2 

P - P - reported  cell  depression  pressure  in 
S1  s2  inches  Hg  or  inches  H2 0. 


o subscript  used  to  denote  ambient  conditions. 

R subscript  used  to  denote  area  around  the  engine 
air  intake. 

Because  V = 0 and  the  energy  change  due  to  changes  in 
elevation  are°neglibible,  the  equation  can  be  rewritten 
in  the  form: 

p - p _ '* 

S1  S2  VR 

0 = 2g^“  losses  (3) 

Losses  are  the  result  of  friction  at  the  surfaces, 
and  turbulence  created  at  expansions,  contractions,  and 
turns . 
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To  evaluate  the  losses,  It  is  necessary  to  have  a good 
physical  description  of  the  cell.  It  must  include  the 
lengths  shown  as  L..  and  L?/  the  cross-sectional  area  of 
the  inlet  stack,  tne  shape  and  dimension  of  any  noise 
suppression  baffles  in  the  inlet  stack,  and  the  cell 
cross-sectional  area  at  the  engine  inlet.  These  can  be 
found  in  engineering  drawings  of  the  test  cell,  or  by 
measuring  the  cell.  Appendix  B includes  a sample  test 
cell  evaluation  sheet. 

The  loss  terms  that  appear  significant  are  the  entry 
losses,  losses  due  to  flow  through  the  baffles,  expansion 
from  the  baffles  into  the  entry  stack,  and  the  flow  through 
the  elbow.  The  friction  losses  through  any  unobstructed 
areas  are  two  to  three  orders  of  magnitudes  smaller  than 
the  above.  All  of  these  losses  can  be  readily  evaluated 
with  the  aid  of  the  Chemical  Engineers'  Handbook  (Ref- 
erence 12)  , 

The  calculation  requires  a measured  cell  depression 
pressure  and  a Reynolds  number  (NR  ) based  on  estimated 
velocity.  Cell  depression  is  readily  obtained  at  various 
fuel  flows  for  various  engines  from  cell  operational  data. 
The  cell  depression  used  in  the  equation  must  be  an  actual 
cell  depression  for  the  specific  test  cell  fuel  flow,  and 
engine  combination  being  evaluated.  An  initial  velocity 
of  20  feet  per  second  is  a good  approximation  for  the  N 
calculation.  A second  iteration  of  the  equation  may  beKe 
used  if  20  ft/sec  is  too  far  off. 

Once  the  reference  velocity  VR  is  known,  it  can  be 
multiplied  by  the  reference  cross-sectional  area  to  get 
the  volumetric  flow  in  ft3/sec.  Using  the  density  of  air 
at  the  cell  conditions  (Reference  11) , the  molecular  weight 
of  air,  and  the  ideal  gas  law,  the  molar  flow  rate  can  be 
determined.  This  number,  plus  the  number  of  moles  of 
gaseous  products  produced  by  the  engine  combustion  yields 
the  molar  flow  rate  at  the  exit  plane  of  the  exhaust  stack. 
Using  the  measured  exhaust  plane  temperature,  atmospheric 
pressure,  and  the  ideal  gas  law,  the  volumetric  flow  rate 
of  the  gases  in  cubic  meters  per  second  can  be  determined. 
The  increase  in  total  flow  rate  due  to  engine  combustion 
gases  is  about  three  orders  of  magnitude  less  than  total 
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cell  air  flow  rate.  The  air  flow  rate  calculated  above 
may  be  used  directly  as  total  cell  mass  flow  rate.  Divid- 
ing the  carbon  production  rate  by  the  test  cell  air  flow 
rate  results  in  a soot  loading  in  g/m3  at  the  exit  plane  of 
the  cell,  at  actual  conditions. 

4.  DETERMINATION  OF  THE  PLUME  OPACITY 


The  estimation  of  plume  opacity  is  made  by  using  the 
Beer-Lambert  Law 

- (S  uiL) 
t = exp  ' p 


where 

t is  the  transmittance  through  the  plume 
Sp  is  the  specific  projected  particle  extinction 
a)  is  the  soot  loading  in  g/m3  at  the  stack  exit 
L is  the  path  length  of  the  beam. 


The  path  length  is  obtained  from  the  geometry  of  the 
cell,  and  it  is  the  longer  of  the  two  stack  dimensions. 

The  longer  dimension  must  be  used  so  that  the  worst  possible 
case  is  considered. 


The  selection  of  a value  for  S is  a difficult  problem. 
McDonald  (Reference  13)  reports  paFticle  diameters  of  O.Olu 
to  0.05u  leaving  the  engine.  Stockham  and  Betz  (Reference 
6)  report  that  the  measured  geometric  mean  particle  diameter 
is  0.05y«  Using  O.Olu  and  0.05u  as  the  two  extremes,  and 
assuming  monodispersed  aerosols,  a Brownian  agglomeration 
calculation  was  done  (Appendix  A) . This  shows  that  agglomer 
ation  could  reduce  the  number  of  particles  by  a factor  of 
up  to  two.  Agglomeration  fox  the  same  mass  concentration 
of  polydispersed  particles  in  the  same  size  range  would  be 
considerably  greater.  Because  of  these  factors,  it  is 
extremely  difficult  to  predict  particle  size  at  the  test 
cell  exhaust. 


Actual  data  for  an  evaluation  of  S is  extremely 
limited.  The  only  data  known  to  this  author  is  in  work 
reported  by  Grems  (Reference  9) . His  work  measuring  the 
particle  size  distribution  at  the  exhaust  plane  of  a test 
cell  showed  a bimodal  particle  distribution  with  a sub- 
micron mode  in  the  range  of  0.08u  to  0.2u  and  then 


12 


larger  particles  in  the  5u  and  greater  range.  The  larger 
particles  are  thought  to  be  the  result  of  scouring  within 
the  test  cell. 

Some  actual  measurements  of  opacity  versus  soot  load- 
ing for  J-57  engines  are  also  reported  by  Grems . Using 
known  values  for  t , ui,  and  L,  values  for  S were  calcu- 
lated. For  low  smoke  numbers  (SN  <25 ) , S ” was  in  the 
range  of  1 to  4m2/g  and  for  SN  >25,  S was  in  the  range 
of  7 to  9 m2/g.  Based  on  this  information  the  following 
values  for  are  suggested: 

SN  <25  S = 3 m2/g 

IT 

SN  >25  Sp  = 8 m2 /g 

These  values  are  reasonable, based  on  the  limited  actual 
data.  Predicting  values  using  the  k parameter  of  Ensor 
and  Pilot  requires  a knowledge  of  geometric  mean  radius 
and  geometric  standard  deviation,  and  actual  physical 
characteristics  of  the  agglomerated  particles.  Because 
these  cannot  be  readily  predicted,  the  use  of  the  empirical 
results  is  recommended. 
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SECTION  III 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  It  is  possible  to  predict  plume  opacity  from  known 
engine,  fuel,  and  test  cell  parameters. 

2.  For  the  data  set  used,  good  correlation  is  shown 
between  predicted  and  actual  soot  loadings.  However, 
problems  could  occur  because  of  the  difficulty  in  pre- 
dicting the  specific  projected  particle  extinction  area 
for  particles  at  the  test  cell  exhaust  stack. 

3.  Figures  to  evaluate  the  test  cell  predicted  flow 
rates  with  actual  flow  rates  were  not  available.  Actual 
cell  flow  rates  should  be  compared  with  the  calculated 
values  when  data  are  available. 

4.  The  value  assigned  to  Sp  for  the  final  calculation  of 
transmittance  was  the  most  likely  source  of  error.  Actual 
soot  loadings  and  opacity  readings  were  used  to  estimate 
the  value  for  S , but  there  were  so  few  actual  readings 
that  the  data  base  was  not  considered  adequate.  It  is 
anticipated  that  a data  base  will  become  available  in 

the  foreseeable  future,  and  when  it  does,  S must  be 
reevaluated.  p 

5.  The  estimation  of  plume  opacity  as  described  in  this 
work  will  provide  a reasonable  first  evaluation  of  test 
cell/engine  combinations  that  could  exceed  visible  emission 
standards.  It  can  also  be  used  to  develop  priorities  for 
corrective  action,  and  establish  the  combinations  of  en- 
gines and  cells  which  would  exceed  standards.  By  shifting 
certain  engines  to  different  cells,  violations  of  emission 
standards  might  be  reduced  with  no  controls  required. 
Utilization  of  this  type  of  evaluation  should  minimize  the 
cost  of  bringing  USAF  Test  Cells  into  compliance  with 
visible  emission  standards. 


SECTION  IV 


SAMPLE  CALCULATION 


Calculation  of  Soot  Loading  at  Engine  Exhaust  Plane. 

Parameters 


Engine : 

Uncorrected  Barometric  Pressure: 
Inlet  Air  Temperature: 

Fuel  Flow 

Smoke  Number  (SAE  SN)  ; 

Weight  per  unit  Area  (W/A) : 

Fuel  to  Air  Ratio: 

1.  Soot  Loading  (SN  > 15) 


J-57-P21 
28.72  in  Hg 
52  °F 

5340  lb/hr 
59.5 

0.0230  lb/in2 
(1.62  g/cm2) 

0.0117 


UJ 

CO 

(0 


e/»e  x W/A  = lo10' 026  3 (S1))  ~ 5 ' 961 g/cm1 

e/»e  - l0«>. 0263  (59. 5l-5.96>g/om>  x onVl.62  g 


n — 5 

e/(e  = 2.49  x 10  parts  soot/parts  exhaust  (by  wt) 


2 . Engine  Flow  Rates 

5340  lb/hr  x 454  g/lb  x .1  hr/3600  sec  x 1 mole/120  g s 
5.61  moles/sec  (fuel  flow) 

5.61  moles  fuel/sec  x 4.7  moles  products/mole  fuel  - 
26.4  moles  products/sec  (product  gas  flow) 
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0.0117  lb  fuel/lb  air  = 5340  lb  fuel/x 
x = 4.56  x 105  lb  air/hr 

4.56  x 105  lb  air/hr  x 1 hr/3600  sec  x 454  g/lb  x 

1 mole/28.97  g = 1985  moles  air/sec  (engine  air  flow) 
Engine  Mass  Flow 

1985  moles  air/sec  + 26.4  moles  products/sec  = 

2011  moles/sec 

2011  moles  exhaust/sec  x 28.95  g/mole  = 

5.82  x lO^g  exhaust/sec 

Exhaust  particulate  flow 

5.82  x 10‘*g  exhaust/sec  x 2.49  x 10  5g/mole  = 

1.45  g soot/sec 

Assuming  no  losses  in  the  test  cell,  soot  loading  at 
the  exit  plane  of  test  cell  exhaust  will  be  1.23  g/sec. 

3.  Air  Flow  Through  the  Test  Cell 

Figure  3 represents  a typical  test  cell. 

(English  units  are  used  because  current  engineering  hand- 
books are  in  English  units.) 

Sample  Dimensions 

Inlet  Area  = 35  ft  x 30  ft  = 1050  ft2  (overall  inlet  area) 

Effective  Area  = Ag  = 700  ft2  (Reduced  by  baffles) 
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31  baffles,  4 in  wide,  35  ft  long  mounted  on 
1 ft  centers. 

Creating  30  parallel  channels,  8 in  wide,  35  ft  long, 

15  ft  deep.  (L^) 

Area  E = 30  ft  x 35  ft  = 1050  ft2  = 1050  ft2 (elbow 
entrance) 

(Flow  in  this  area  will  be  in  a square  or  mitered  elbow.) 
Area  R is  30  ft  x 30  ft  = 900  ft2 . 

Cell  Depression  Pressure  = 0.25  in 
Velocity  Relationships 

X AB  * VE  X AE  = X AR 

VB  s VR  Ar/Ab  = (900  ft2  /700  ft2)  x VR  - 1.29  VR 
VE  = VRAR/AE  = (900  ft2  /10 50  ft2)  x VR  = . 857  VR 

Losses : 

a.  Entry 

Fent  = Kent  x V/2gc*  Since  VAo  s 700  ft2/°°  = 0 ' 

K = 0.5  for  turbulent  flow  from  Chem  Engineers 

ent  ,,  ,,  . 3 

Handbook 

Fent  = (0.5)  ( 1 . 29 ) 2 x VR2/2gc  = 832  VR2/2gc 

b.  Flow  Through  Baffles 

Hydraulic  Radius  = Area  of  Stream  Cross-Section/Wetted 

Perimeter  (for  1 channel) 

Hydraulic  Diameter  (DR)  = 4^ 
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Dimension  of  openings  for  4 in  baffles  on  12  in  centers 
= 8 in  wide  by  35  ft  long 

Area  of  Cross-section  = 0.667  ft  x 35  ft  = 23.3  ft2. 

Wetted  Perimeter  = 2 (0.667  ft)  + 2 (35  ft)  = 71.3  ft. 

Rh  = 23.3  ft2 /71. 3 ft  = 0.327  ft.  DR  = 4RH  = 1.31  ft. 

NRe  * Reynolds  number  = DR  V p/p (Approximate  Velocity 
is  20  ft/sec) . 

p = air  density  = 0.0687  lb/ft3  at  28.72  in  Hg  and 

52  0 F.  (9) 

- 5 

p = air  viscosity  = 1.21  x 10  lb/'(ft)(sec)  (Chem  Eng 
Handbook ) . 

NRe  s (1.31  ft)  (20  ft/sec)  x 0.0687  lb/ft3  x (ft) (sec)/(1.21 
x 10“5)  lb  * 1.5  x 105. 

e = surface  roughness  = 0.01  ft  (Rough  concrete  is  a 

good  approximation  of  the  roughness  of  the  baffles. 
However,  it  may  be  necessary  to  calculate  a better 
value  if  the  baffles  are  extremely  rough) . 

e/DH  ■ 0,01  ft/1.31  ft  = 0.008 

Use  e/Dh  (0.008)  and  NRg  (1.6  x 10 6)  obtain  friction 
factor,  f - 0.009  (Chem  Eng  Handbook). 

L = 15  ft  (L-l) 

F - 4 (f)(L)/DH  x V02/2gc  - ( (4)  (0.009)  (15  ft)/{1.31  ft)) 

X (1.29) 2 VR2/2gc 
F * 0.686  VR2/2gc 
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c. 


Expansion  from  Baffles  to  Elbow 
F - y/2gc  X U-ab/ae>? 

F = (1.29)2  VR  2/2gc  x (1-700/1050  )2 

F = 0.185  VftV2gc 

d.  Flow  Through  Elbow 

F = x VE2/2gn  Kel^  =1.3  (Chem  Eng  Handbook) 

F = (1.  3)  (0.857)  2 V 2 /2g 

K C 

F = 0.955  VR2/2gc 

e.  Flow  to  Engine 

30  ft  x 30  ft;  cross-section  AR  =•  900  ft2 

Wetted  Perimeter  = (4) (30  ft)  = 120  ft 

Rh  = 900  ft2/120  ft  = 7.5  ft;  DR  = (4) (7.5  ft) 

* 30  ft 

Use  same  p,  V,  and  p as  for  baffles. 

NRe  = (30  ft) (20  ft/sec)  x (0.0687  lb/ft3) 
x (ft  sec/1.21  x 10“5lb)  - 3.41  x 106 
Use  t for  smooth  concrete  0.001  ft 

e/D  = 0.001  ft/30  ft  = 3.33  x lO-5;  L = 13  ft(L,,) 
f s 0.0025  (from  chart,  Chem  Eng  Handbook). 

F - (4)  (0.0025)  (13  f t)  / ( ( 30  ft)  x (VR2  /2gc)  ) 

F = 0.0043  VRz/2gc  (negligible  compared  to 
other  losses). 
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Finding  VR 

From  Eq  (3),  P - P = Reported  Cell  Depression 
S1  s2 

Pressure  = 0.25  in  H20 

0 - PSl  - V»  ' (V/29c>  - (»-832VRV2gc) 

- (0.686  V„!/2g„)  - (0.185  v’/2gj 

K C K C 

- (0.955  VR2/2gc) 

0 = (0.25  in  H20  ft3/0.0764  lb) (5.2  lb/£t2in  H20) 

- (3.66  VR2/<2  x 32.17  ft  sec2)) 

0 = 17.0  ft  - (0.0569  ft/sec2 ) V 2 

K 

VR  = 299  ft2 /sec2 

VR  = 17  ft/sec.  Since  20  ft/sec  was  assumed  for  VR, 

VB  = 1.29  VR  = 21  ft/sec 

Second  calculation  not  required. 

Calculation  of  Volume  Flow  Rate  i .. . cell 

17  ft/sec  x 900  ft2  x (lm/3.28  ft) 3 - 430  m3/sec 

(additional  volume  of  combustion  products  is 
negligible  by  comparison.) 
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Volume  flow  at  exit  stack 

Inlet  temperature  ■=  52°F  - 284°K 

Exit  temperature  = 120 °C  = 303°K 

Inlet  pressure  = 28.70  (barometric  pressure  minus 

cell  depression) 

Exit  pressure  = 28.72 

430  m3/sec  x (393/284)  x (28.70/28.72)  = 595  m3/sec 
Soot  Loading 

w = (1.45  g/sec)  x (sec/595  m3)  = 2.44  x 10~3 g/m3 
Plume  Opacity 

x = exp  - ((8  mz/g)  x (2.44  x 10  3g/m3)  x (10.7  m) 
t = 0.81;  Opacity  = 19  percent. 
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APPENDIX  A 


AGGLOMERATION  CALCULATION  (13) 

Length  of  diffuser  - approximately  30  feet 

Exhaust  gas  exit  velocity  at  5000  lb/hr  fuel  flow  = 

1500  ft/sec 

- 2 3 

Soot  Loading  2.601  x 10  g/m 

-6  -4 

Particle  diameter  = 0.02w  or  0.02  x 10  m or  0.02  x 10  cm 
Exhaust  gas  temperature  = 900 °F 

3 

Density  of  carbon  = 2 g/cm 
Volume  of  a particles 

V = 4/3  r3 

V = 4/3  x it  x (0.01  x 10*"*cm)  3 
v - 4.19  x 10  18cm3/particle 

Particles  departing  the  engine: 

Nq  = (particle/4.19  x 10~38cm3)  x (cm3/2g)  x 

(2.601  x 10-2g/m3)  x (m3/ ( 102cm> 3) 

9 3 

Nq  * 3.10  x 10  particles/cm 

Evaluation  of  K0,  coagulation  constant 

Fuchs  gives  KQ  for  a 0.01  radius  particle  to  be 

34  x 10~^°cm3/sec.  It  is  computed  at  a temperature  of 
60 °F  or  520 °R.  To  adjust  for  the  exhaust  gas  temperature  of 
1360  °R, 

Ta/To  . Ko  e 1360°R/520°R  x 34  x 10"10cm3/sec  « 


1 


Evaluation  of  time  in  diffuser 


The  velocity  at  the  exit  of  the  engine  is  1500  ft/3ec. 

The  gas  is  slowed  through  the  diffuser  so  that  a velocity 
of  1000  ft/sec  will  give  a good  estimate  of  residence  time 
of  the  gas  in  the  diffuser. 

30  ft  x sec/1000  ft  = 0.03  second  in  diffuser 

Evaluation  of  agglomeration: 

From  Fuchs: 

N = N /(I  + K N t) 
o o o 

N = (3.10  x 10 9 particles/cm3) / (1  + ((8.89  x 10~9cm3/sec)  x 
9 3 

(3.10  x 10  particles/cm  ) x 0.03  sec) 

N = (3.10  x 109  particles,/cm3)  / (1  + 0.827) 

N ~ 1.70  x 109  particles/cm3 

N/No  = (1.70  x 109  particles/cm3) /( 3. 10  x 109  particles/cm3) 

= 0.55.  Thus  the  number  of  particles  is  reduced  by 
approximately  a factor  of  two. 

Doing  the  same  calculation  using  a particle  diameter  of  0.05 
Volume  of  particle: 

V = 4/3  x v x (0.025  x 10-4cm)3 

V = 6.54  x 10  ^cm3/particle 
Particles  departing  engine: 

Nq  = (particle/6.54  x lo"17cm3)  x (cm3/2g)  x (2.601  x 
10'2g)  x (m3/ (102cm) 3) 

Nq  = 1.98  x 10®  particles/cm3 
Evaluation  of  K 

o 

Kq  is  a function  of  temperature:  (the  same  as  for  .01  ) 

Kq  » 8.89  x 10~9cm3/sec 
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Evaluation  of  time  in  diffuser: 

t will  not  change  unless  fuel  flow  changes 
Evaluation  of  agglomeration: 

N ■=  (1.98  x 10®  particles/cm®) /( 1 + (8.89  x 10  ^cnw'sec)  x 
8 3 

(1.98  x 10  particles/cm  ) x 0.03  sec) 

N - 1.88  x 10®  particles/cm® 

N/No  - 0.95  (no  significant  agglomeration) 


APPENDIX  B 


TEST  CELL  INFORMATION  SHEET 


1.  The  information  requested  will  be  used  to  assess  the 
environmental  impact  of  test  cell  operations  at  your  base. 

The  information  is  required  only  once  for  each  type  of  test 
cell . 

2.  You  may  use  schematic  drawing  (Figure  B-l)  of  a typical 
test  cell  to  supply  required  test  cell  dimensions.  Simply 
write  in  the  dimension  indicated  beside  each  of  the  double 
headed  arrows.  Cross-section  A-A  on  the  schematic  corres- 
ponds to  the  overall  area  of  the  inlet  stack  entrance. 
Cross-section  B-B  is  the  area  of  the  cell  exhaust  stack. 
Cross-section  C-C  corresponds  to  the  area  around  the  cell 
depression  pressure  probe  (usually  at  the  engine  intake) . 

3.  If  the  cell  has  any  special  features  or  if  it  differs 
from  the  sketch,  feel  free  to  modify  the  sketch  as  required, 
labeling  the  changes  clearly. 

4.  In  addition  to  the  dimensions  requested  on  the  sketch, 
please  answer  the  following  questions: 

a.  How  is  the  cell  identified?  (e.g..  Bldg  Number, 
site  number,  etc.) 

b.  What  type  of  engines  are  run  in  this  cell?  List 
each  type.  (e.g.,  J-57-P21  A) 

c.  How  far  in  front  of  the  engine  is  the  cell  depres- 
sion probe?  Indicate  on  schematic  diagram. 

d.  Where  is  the  compressor  inlet  temperature  probe? 
Indicate  on  schematic  diagram. 

e.  How  many  cooling  water  rings  are  there?  Where  are 
they  located?  (Mark  on  sketch  W..,  W«,  W^,  etc.)  and  give 
distances  from  an  identified  point. 

f.  Does  the  air  intake  have  baffles?  If  possible, 
send  a copy  of  the  plans  showing  the  baffles. 

g.  If  the  baffle  plans  are  not  available,  please  make 
a sketch  showing  their  locations  in  the  intake,  the  shape  of 
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the  baffles,  and  their  dimensions.  Be  sure  to  include  the 
followings 


(1)  Number  of  openings  for  inlet  air  (spaces 
between  baffles) . 

(2)  Length  and  width  of  each  baffle  opening. 

(3)  Depth  (how  far  do  baffles  extend  into  inlet 

stack?) 

h.  How  far  below  the  exhaust  stack  exit  is  the  exhaust 
gag  temperature  probe? 

i.  Line  A-B  on  the  sketch  represents  a small  wall  that 
sometimes  protrudes  into  the  horizontal  portion  of  the  cell 
from  the  top,  sides,  and/or  floor.  This  is  normally  between 
the  intake  stack  and  horizontal  portion  of  the  cell.  If 
there  is  a restriction  of  this  type,  give  us  the  cross- 
sectional  dimensions  of  the  cell  at  this  point,  as  well  as 
at  the  cell  depression  probe. 

i Restricted  Height  Width 

I , 

I 1 

| Thickness  of  Protrusion  

Again,  the  cell  information  requested  must  only  be 
provided  once.  Please  try  to  make  it  as  accurate  and  complete 
as  possible.  In  particular,  the  cell  depression  pressure 
must  be  measured  as  carefully  as  possible.  The  decision  to 
implement  a smoke  reduction  program  for  your  test  cells 
depends  largely  on  this  information. 
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APPENDIX  C 


FLOW  CHART  FOR  CALCULATION  OF  PREDICTED  OPACITY 


The  simplified  flow  chart  on  the  following  pages  indicates 
the  most  efficient  method  of  computing  the  predicted  opacity 
for  a given  test  cell  and  engine  combination.  This  general 
scheme  was  used  in  the  development  of  a computer  program 
written  in  FOCAL  for  the  PDP-8  minicomputer.1 

The  computations  are  divided  into  two  main  parts.  Since 
the  friction  losses  for  a given  test  cell  are  constant  regard- 
less of  the  engine  used,  these  are  calculated  first.  Using 
test  cell  dimensions,  the  loss  is  computed  for  each  contraction, 
expansion,  baffle,  and  elbow  encountered  by  the  air  as  it  flows 
from  the  test  cell  inlet  to  the  engine  intake. 

The  total  friction  losses  are  then  used  to  calculate  the 
velocity  and,  hence,  volume  of  air  flow  in  the  reference  section 
of  the  *-est  cell,  using  the  measured  cell  depression  for  a given 
engine  operation.  The  other  measured  parameters  for  that  engine 
operation,  i.e.,  SAE  SN , Fuel  Flow,  Fuel/air  ratio,  are  used  to 
calculate  the  particulate  concentration  in  the  test  cell  exhaust 
which,  in  turn,  is  used  to  compute  a predicted  opacity. 

The  detailed  calculations  referred  to  .in  the  flow  chart 
may  be  found  in  Section  IV  of  this  report,  "Sample  Calculations." 


* A copy  of  the  program  may  be  obtained  by  contacting  the 
Division  Chief,  AFCEC/EVA,  Tyndall  AFB  FL  32403. 
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MAIN  FLOW  CHANT 
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Elbow  Friction  Loss  Procedure 


l 


Expansion  Friction  Loss  Procedure 


3 3 


Reference  Velocity  Computation 


Use  total  friction  loss,  air  density, 
cell  depression  pressure  and  reference 
section  cross-sectional  area  to  compute 
volume  flow  in  reference  section  in 

m3/sec , 


Exit  Stack  Volume 


Flow 


Use  compressor  inlet  temperature,  exit 
stack  temperature,  exit  stack  pressure, 
and  cell  depression  pressure  to  compute 
volume  flow  at  exit  stack. 


Particulate  Mass  Flow 


Use  SAE  SN,  fuel  flow,  and  fuel/air  ratio 
to  compute  particulate  mass  flow  at  engine 
exhaust  in  g/sec. 


Predicted  Transmittance 


Use  stack  exit  volume  flow  and  particulate 
mass  flow  to  determine  particulate  mass 
concentration  in  g/m3.  Use  this  value  along 
with  plume  width  to  calculate  percent  trans- 
mittance . 
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INITIAL  DISTRIBUTION 


HQ  USAF/PREV-P  1 
HQ  USAF/PREV-X  1 
HQ  USAF/RDPS  2 
HQ  USAF/SAFOI  1 
CINCAD/DEEV  2 
AFLC/SGB  1 
AFLC/DEPV  1 
AFLC/MAUT  1 
AFLC/MMRF  1 
AFSC/DE  1 
AFSC/SD  1 
AF SC/DEV  1 
AFSC/SGB  1 
ATC/DEPV  1 
ATC/SGPAP  1 
MAC/SGPE  1 
MAC/DEEE  1 
MAC/DEMP  1 
CINCSAC/DEPA  2 
CINCSAC/DEPV  1 
CINCSAC/SGPA  1 
TAC/DE  1 
TAC/DEEV  1 
TAC/SGPB  1 
AFRES/DEEE  I 
AUL  (AU-LSE- 70-239)  1 
AFOSR  1 
AFAPL  1 
AFML/DO  (Library)  1 
OEHL/CC  3 
OEHL/OL-AA  1 
OEHL/OL-AB  1 
AFWL/SUL  (Tech  Library)  1 
AFRPL/ Library  1 
FTD/LGM  1 
AEDC/DEE  1 
ADTC/CSV  1 
ADTC/DLOSL  (Tech  Library)  1 
USAF  Rgn  Civ  Engrg/Atlanta  1 
USAF  Rgn  Civ  Engrg/San  Francisco  1 
USAF  Rgn  Civ  Engrg/Dallas  1 
SAM4A/MAGCB  1 
DDC/TCA  12 
ARPA  1 
Defense  Research  6 Engrg/AD  (EgLS)  1 
OASD/(If,L)  ES  1 
USA  CERL  1 


Dir,  USA  Eng  R5D  Lab/MERDC  1 

Ch  of  R(jD,  Dept  of  the  Amy/ 
DARD-ARE-E  1 

AFRCE/Eastem  Region  1 

AFRCE/Central  Region  1 

AFP.CE/Westem  Region  1 

Chief  of  Naval  Op 

Environ  Protection  Div,  OP- 4 5 1 

NCEL,  Code  2S111  1 

Naval  Air  Dev  Ctr/MAE  1 

Technology  Transfer  Staff  (EPA)  1 
USAFA/DFCBS  20 

National  Science  Foundation  1 

US  Army  Med  Bioengrg  RfjD  Lab 
ATTN:  SGRD-UBG  1 

Air  01B,  NAVAIR  1 

Naval  Postgraduate  School  1 

NARD  (Code  64270)  1 

NAVFAC  ENGCCM  1 

Naval  Weapons  Center  1 

University  of  New  Maxico  1 

AFCEC/SU  1 

AFCEC/SUL  1 

AFCEC/WE  1 

AFCEC/EV  8 
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